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Abstract 

In this technical report, we consider a nonlinear 4th-order degenerate parabolic 
partial differential equation that arises in modelling the dynamics of an incompressible 
thin liquid film on the outer surface of a rotating horizontal cylinder in the presence 
of gravity. The parameters involved determine a rich variety of qualitatively different 
flows. Depending on the initial data and the parameter values, we prove the existence 
of nonnegative periodic weak solutions. In addition, we prove that these solutions 
and their gradients cannot grow any faster than linearly in time; there cannot be a 
finite-time blow-up. Finally, we present numerical simulations of solutions. 
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1 Introduction 

We consider the dynamics of a viscous incompressible fluid on the outer surface of a horizontal 
circular cylinder that is rotating around its axis in the presence of gravity, see Figure [T| If 
the cyfinder is fully coated there is only one free boundary: where the fiquid meets the 
surrounding air. Otherwise, there is also a free boundary (or contact fine) wfiere the air and 
fiquid meet the cyfinder's surface. 

The motion of the liquid film is governed by four pfiysical effects: viscosity, gravity, surface 
tension, and centrifugal forces. Tfiese are refiected in tfie parameters: R — tfie radius of tfie 
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Figure 1: Liquid film on the outer surface of a rotating horizontal cylinder in the presence 
of gravity. 

cyhnder, u — its rate of rotation (assumed constant), g — the acceleration due to gravity, u 
— the kinematic viscosity, p — the fluid's density, and a — the surface tension. 

These parameters yield three independent dimensionless numbers: the Reynolds number 
Re = {R^Lo)/v, the Galileo number Ga = g/{Ruj'^) and the Weber number We = {pR^uj'^)/a. 

We introduce the parameter e = h/R, where h is the average thickness of the liquid. The 
following quantities are assumed to have finite, nonzero limits in the thin film (e — > 0) limit 
[271 EH El [23]: 

K = Ree^, X = =^e^5 and /i = GaRee^. (1.1) 

This corresponds to a low rotation rate, for example. 

One can model the flow using the full three-dimensional Navier-Stokes equations with 
free boundaries: for u{x,y,z,t) in the region x G [— vr, tt), y G M^, and z G {0, h{x,y,t)) 
where x is the angular variable, y is the axial variable, and h{x,y,t) is the thickness of the 
fluid above the point (x, y) on the surface of the cylinder at time t. This has been done 
by Pukhnachov |27| in which he considered the physical regime for which the ratio of the 
free-fall acceleration and the centripetal acceleration is smaU. There, he proved the existence 
and uniqueness of fully-coating steady states (no contact line is present). We know of no 
results for the affiliated initial value problem. 

In this physical regime, if one also makes a longwave approximation (the thickness of 
the coating fluid is smaUer than the radius of the cylinder) and if one further assumes that 
the rotation rate is low (or the viscosity is large) then the three-dimensional Navier-Stokes 
equations with free boundary can be approximated by a fourth-order degenerate partial 
differential equation (PDE) for the fllm thickness h{x,y,t). This is done by averaging the 
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fluid flow in the direction normal to the cylinder (23 [2B]- If one further assumes that the 
flow is independent of the axial variable, then this results in a PDE in one dimension for 
h{x, t). 

In his pioneering 1977 article about syrup rings on a rotating roller, Moffatt neglected 
the effect of surface tension (i.e. We~^ = = x)? assumed the flow was uniform in the axial 
variable, and derived |23| the following model for the thin film thickness: 

(/i- f/i^ cos(x))^ = 0, (1.2) 



where is given in (1.1) and 

X G [— 7r,7r], t > 0, h is 27r-periodic in x. 
Pukhnachov's 1977 article [57] gives the first model that takes into account surface tension: 

ht + ih- fh^ cos(x)),. + f {h^ (K + K..)), = (1.3) 



where /i and x are given in (1.1) and 

X G [— 7r,7r], t > 0, his 27r-periodic in x. 
This model assumes a no-shp boundary condition at the hquid/solid interface. For a solution 



to (1.2) or (1.3) to be physically relevant, either h is strictly positive (the cylinder is fully 



coated) or h is nonnegative (the cylinder is wet in some region and dry in others). 



Surprisingly little is understood about the initial value problem for (1.3). Bernis and 
Friedman [3] were the first to prove the existence of nonnegative weak solutions for nonneg- 
ative initial data for the related fourth-order nonlinear degenerate parabolic PDE 

ht + U\^)h...). = Q, (1.4) 

where /(/i) = |/ir/o(/i), /o(/^) > 0, n^l. 

Unlike for second-order parabolic equations, there is no comparison principle for equation 



(1.4). Nonnegative initial data does not automatically yield a nonnegative solution; indeed 



it may not even be true for general fourth-order PDE (e.g. consider ht = —hxxxx)- The 



degeneracy f{h) in equation (1.4) is key in ensuring that nonnegative solutions exist. Also, 
unhke for second-order problems, it is possible that strictly positive initial data might yield 
a solution that is zero at certain moments in time, at certain locations in space. 



3 



Lower-order terms can be added to equation (1.4) to model additional physical effects. 
For example, 

ht + {f{h)h^xx)x-{g{h)h^)^ = (1.5) 



where g{h) > for /i 7^ 0. Equation (1.5) can model a thin liquid film on a horizontal surface 
with gravity acting towards the surface. If this surface is not horizontal then the dynamics 
can be modelled by 



ht + {K^{a -bh^ + h^^^))^ = 0, a > 0, 6 > 



(1.6) 



The constant a in the first-order term vanishes as the surface becomes more and more 
horizontal. If the thin film of liquid is on a horizontal surface with gravity acting away from 
the surface then the thin film dynamics can be modelled by 



ht + {f{h) h^^^)^ + {g{h)h,j,)^ = 0. 
For a thorough review of the modelling of thin hquid films, see |13l EH 126] . 



(1.7) 



In equations (1.5) and (1.6) the second-order term is stabilizing: if one linearizes the 



equation about a constant, positive steady state then the presence of the second-order term 



increases how quickly perturbations decay in time. In equation (1.7), the second-order term 



is destabilizing: the linearized equation can have some long-wavelength perturbations that 



grow in time. For this reason, we refer to equation (1.7) as "long-wave unstable". The long 



wave stable equations (1.5) and (1.6) have similar dynamics as equation (1.4) however the 



long- wave unstable equation (1.7) can have nontrivial exact solutions and can have finite- 



time blow-up {h{x*,t) I 00 as t t t* < 00). 

In all cases, the fourth-order term makes it harder to prove desirable properties such as: 
the short-time (or long-time) existence of nonnegative solutions given nonnegative initial 
data, compactly supported initial data yielding compactly supported solutions (finite speed 
of propagation), and uniqueness. Indeed, there are counterexamples to uniqueness of weak 



solutions [3]. Results about existence and long-time behavior for solutions of (1.5) can be 



found in |5j; analogous results for (1.6) are in [T^. See [HI E] for results about existence. 



finite speed of propagation, and finite-time blow-up for equation (1.7) 



In this paper we study the existence of weak solutions of the thin film equation 

ht + (|/i|^(ao h^xx + a-ihx + a2w'{x)))^ + a^^h^ = (1.8) 
where Oi, 02, 03 are arbitrary constants, constant ao > 0, and w{x) is periodic. Equation 



(1.3) is a special case of (1.8). The sign of ai determines whether equation (1.8) is long-wave 
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unstable. Also, the coefficient of the convection term a2{w'{x)\hf)x can depend on space and 



will change sign if a2w'(a;) ^ 0. The cubic nonlinearity \h\ in equation (1.8) arises naturally 



in models of thin liquid films with no-shp boundary conditions at the hquid/solid interface. 
Our methods generalize naturally to f{h) = for an interval of n containing 3; we refer 
the reader to [31 El E] for the types of results expected. 

Given nonnegative initial data that satisfies some reasonable conditions, we prove long- 
time existence of nonnegative periodic generahzed weak solutions to the initial value problem 



for equation (1.8). We start by using energy methods to prove short-time existence of a weak 



solution and find an exphcit lower bound on the time of existence. A generalization and 
sharpening of the method used in [S] allows us to prove that the norm of the constructed 
solution can grow at most hnearly in time, precluding the possibility of a finite-time blow-up. 
This control, combined with the explicit lower bound on the (short) time of existence, 
allows us to continue the weak solution in time, extending the short-time result to a long-time 
result. 

If 02 = or 03 = in equation (1.8) then solutions will be uniformly bounded for all 



time. If 02 7^ and 03 7^ 0, it is natural to ask if the nonlinear advection term could cause 
finite-time blow-up (h{x*,t) | 00 as t j t* at some point x*). Such finite-time blow-up is 
impossible by the linear-in-time bound on but we have not ruled out that a solution 
might grow in an unbounded manner as time goes to infinity. 



In [mdl], the authors consider the multidimensional analogue of (1.4) 



+ V - (|/irVA/i) =0, 



(1.9) 



for h{x, t) where x G C with G {2, 3}. Depending on the sign of A\ if = then 
equation 

ht + V- if{h)\/Ah + VA{h)) = g{t, X, h, Wh) (1.10) 



on Q is the multidimensional analogue of equation (1.5) or (1.7). In [I5], the authors consider 
the long-wave stable case with g = and power-law coefficients, f{h) = and A'{h) = 
— \h\"^. In |18j, the author considers the Neumann problem for both the long-wave stable 
and unstable cases with the assumption that f{h) > has power-law-like behavior near 
/i = 0, that |A'(/i)| is dominated by f{h) (specifically |A'(/i)| ^ c?o/(^) for some do), and 
that the source/sink term g{t,x, h) grows no faster than linearly in h. In [32l |33l |35], the 



authors consider the Neumann problem for the long-wave stable case of ( |1.10[ ) with power- 



law coefficients and a larger class of source terms: g{t, x, h) ~ \h\^ with A > 0. In 
the same authors consider the long-wave stable equation with power-law coefficients but with 
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g{h) = a- \/b{h) where b{z) ~ and a E M^: g models advective effects. Tfiey consider tfie 
problem both on and on a bounded domain Q. 



All of these works on (1.9) and (1.10) construct nonnegative weak solutions from nonneg- 
ative initial data and address qualitative questions such as dependence on exponents n and 
m and A, on dimension N, speed of propagation of the support and of perturbations, exact 
asymptotics of the motion of the support, and positivity properties. We note that the works 
[321 61 Eg [301 El] also construct "strong" solutions. 



2 Steady state solutions 

Smooth steady state solutions, h{x,t) = h{x), of (Oj) satisfy 



h - f /i^ cos(a;) + f {h^ {K + Kxx)) = q (2.1) 

where g is a constant of integration that corresponds to the dimensionless mass flux. In the 
zero surface tension case (x = 0), steady states satisfy 

h - f /i^ cos(x) = q. (2.2) 

Such steady states were first studied by Johnson [I9] and Moffatt [23]. Johnson proved 
that there are positive, unique, smooth steady states if and only if the flux is not too large: 
< g < 2/{3y/j2). These steady-states are neutrally stable |25|. 

This critical value of q had been flrst observed numerically by Moffatt as a threshold 



between continuous and discontinuous (shock) steady states. Evaluating (2.2) at x = ±7r/2, 
one sees that this limit to the amount of fluid which can be transferred per unit time corre- 
sponds to a limit to the thickness of the fluid at the top (or bottom) of the cylinder. 

Figure [2] presents steady states for two fluxes. The smooth curve corresponds a steady 
state with a flux smaller than 2/{3y/JI). As q decreases to zero, the thickness of the fluid 
on the left side of the cylinder decreases to zero. As q increases to 2/(3^/1) the smooth 
maximum on the right side of the cylinder becomes a corner, as shown in Figure [2j Also, as 
shown, at this critical flux value there can be discontinuous steady states. 

Smooth, positive steady states in the presence of surface tension have been studied by 
a number of authors. One striking computational result [I] is that for certain values of 
X and /i there can be non-uniqueness. Speciflcally, one can flnd flux values q for which 
there are more than one steady state with that flux (the steady states have different total 
mass). Similarly, one can find total masses for which there are more than one steady state 
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Figure 2: Three steady state solutions of equation (2.2) with /i = 1. Left: solutions 
plotted as {x,h{x)). The smooth curve corresponds a steady state with flux q = 0.64. The 
discontinuous curve and the curve with a corner in it correspond to the critical flux value 
q = 2/3. Right: solutions plotted as a film coating a cylinder of radius 1. The cylinder is 
denoted with a heavy fine. 

(the steady states have different flux). These steady states were numerically discovered via 
an elegant combination of asymptotics and a two-parameter (mass and flux) continuation 
method Figure 14]. To start the continuation method, earlier work [2j on the regime in 
which viscous forces dominate gravity was used. There, asymptotics show that for small 
fluxes the steady state is close to g + l/3g^cos(a;) + O(g^), providing a good flrst guess for 
the iteration used to flnd the steady state. The bifurcation diagram shown in Figure 14 of [T] 



also suggests that the Moffatt model (1.2) can be considered as the hmit of the Pukhnachov 
model (1.3) as surface tension goes to zero (x — ^ 0). 

We are not aware of a result that proves that smooth positive steady states exist if and 
only if < g < g*(/i) for some g*(/i). Pukhnachov proved |29] a nonexistence result: no 
positive steady states exist if g > 2^/3//! ~ 3.464/y//i. We improve this, proving that no 
such solution exists if g > 2/3 \/2/ ^ ^ 0.943/y^. 



Pro posit ion 2.1. There does not exist a strictly positive 27r periodic solution h{x) of equa- 



tion (2.1) ifq> 2/3 



Proof of Proposition \2. 1\ FoUowing Pukhnachov, we start by rescaling the flux to 1 by in- 



troducing y{x) = h{x)/q and introducing the parameters 7 = ^ and (3 = Equation 



(2.1 ) transforms to 

7(/' + Z/')=/3cos(x)-^ + 4,. (2.3) 

The solution y is written as 

y{x) = ao + ai cos(x) + 02 sin(a;) + v{x) (2.4) 
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where v{x) ± span{l, cos(x), sin(x)} and satisfies 



(2.5) 



A solution V exists only if the right-hand side of (2.5) is orthogonal to span{l, cos(a;), sin(x)}. 
As a result, 

2n 



1 



1 



dx = 



2lT 



^-^)cos(x)cix = 7r/3, 



(2.6) 



(2.7) 



Adding equations (2.6) and (2.7) yields 

2n 



TT/3 



1 



y{x 



< 



(1 + cos(x)) dx 



The function F{y) = — is bounded above by 4/27 on [1, oo) hence 



Tip < 



27 



[1 + cos(x)) dx < ^271 



This shows that if there is a steady state then (3 < 8/27. Recalling the definition of (3, there 
is no steady state if g > 2/3 \/2/jj,. □ 

The proof also holds in the case of zero surface tension x = 7 = and so it is natural that 
the bound 2/3 \/2/n is larger than 2/(3^/1) (the bound found by Johnson and Moffatt.) 
Also, we note that numerical simulations that suggest nonexistence of a positive steady 
state if g > 0.854 when /x = 1 for a large range of surface tension values |20l P- 61]; 
our bound of 0.943 is not too far off from this. We close the discussion of steady states 
by considering their nonhnear stability. This is done via simulations of the initial value 
problem for different regimes of the PDE. Figure [3] considers the PDE with no advection. 
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ht + {h^{hxxx + 16 /ix))^ = 0. The PDE is translation invariant in x and constant steady 
states are linearly unstable. As a result, any non-constant behaviour observed in a solution 
starting from constant initial data would be due to growth of round-off error. For this 
reason, non-constant initial data is chosen: hQ{x) = 0.3 + 0.02 cos(x) + 0.02 cos(2x). The LF' 
and norms of the resulting solution appear to be converging to limiting values as time 
passes and long-time hmit of the solution appears to be four steady-state droplets of the 
form a cos(4x + 0) + 6 for appropriate values of a, 0, and h. Like the PDE, the simulation 
shown respects the symmetry about x = vr of the initial data. However, we find that if 
one computes longer, the symmetry is broken and the solution appears to converge to a 
profile with three steady droplets. This suggests that the four droplet configuration may be 
a steady state but it's an unstable one and accumulated round-off error eventually leads the 
numerical solution away from it. 

Figure [4] shows the evolution from constant initial data for the PDE with nonlinear 
advection but no hnear advection: ht + {h?{hxxx + IQhx + 9>cos{x)))x = 0. The long- 
time limit appears to be a steady state which is zero (or nearly zercQ on [0, tt] with the 
bulk of the fluid contained in a droplet supported within (tt, 27r), centred roughly about 
the bottom of the cylinder {x = 37r/2). Finally, Figure [s] shows the evolution resulting 
from the same constant initial data for the PDE with both linear and nonlinear advection: 
ht + {h?{hxxx + 16 /ix + 8 cos(x)))x. + S/i^ = 0. The long-time limit appears to be fully wetted 
cyhnder with a steady "droplet" centred shghtly past the bottom of the cylinder (here "past" 
refers to the direction determined by the direction of rotation a;; see Figure [T]). 

We close by noting that the PDE considered in Figure [5] corresponds to coefficient as = 3 



in the PDE (1.8). As we increase the value of 03 we find there appears to be a critical 
value past which the solution appears to converge to a time-periodic behaviour rather than 
a steady state. Specifically, a "thumping" behaviour is observed in which the cylinder is fully 
wetted but the bulk of the fiuid is located in one region. This bulk of fiuid moves around 
the rotating cylinder in a time-periodic manner. 



^The solutions shown have a very thin film of liquid (of order 10 ^) in the apparently dry region. 
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Figure 3: The evolution equation with no hnear or nonlinear advection, ht + {h^{hxxx + 
16/ix))x — 0, corresponding to Oq = 1, Oi = 16, and a2 — — 0. The initial data is 
hQ{x) — 0.3 + 0.02 cos(a;) + 0.02 cos(2x). Left plot: the solution at times t = (dashed line), 
t = 12, 12.5, 13, 15 (solid lines), and t — 140 (heavy line). Right plot: the and norms 
plotted function of time. 




2000 2500 3000 



Figure 4: The evolution equation with nonlinear advection but no hnear advection, ht + 
{h^ihxxx + 16/i3, + 8cos(x)))a: = 0, corresponding to = 1, ai = 16, a2 = 8, and = 0. 
The initial data is ho{x) = 0.3. Left plot: the solution at times t = (dashed line), 
t — 0.5,1,2,10 (sohd lines), and t — 3000 (heavy line). Right plot: the and norms 
plotted function of time. 
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Figure 5: The evolution equation with both hnear and nonlinear advection, ht + {h^ih^xx + 
IQhx + 8cos(a;)))^ + ?>hx = 0, corresponding to = 1, ai = 16, 02 = 8, and 03 = 3. The 
initial data is hQ{x) = 0.3. Left plot: the solution at times t = (dashed line), t = 0.5, 1, 2, 4 
(solid lines), and t = 20 (heavy hne). Right plot: the and norms plotted as a function 
of time. 



3 Short-time Existence and Regularity of Solutions 

We are interested in the existence of nonnegative generalized weak solutions to the following 
initial-boundary value problem: 

ht + {f{h){aohxxx + aihx + a2w'{x)))^ + a^hx = in Qt, (3.1) 
g(-a,t) = g(a,t) for t > 0, ^ = M, (3.2) 
/i(x,0) = /io(a;) ^ 0, (3.3) 

where f{h) = \h\^^ h = h{x,t), VL = {—a, a), and Qt = f2 x (0,T). Note that rather than 
considering the interval (—a, a) with boundary conditions (3.2) one can equally well consider 
the problem on the circle S^; our methods and results would apply here too. Recall that 
ai, a2, and 03 in equation (3.1) are arbitrary constants; oq is required to be positive. The 
function w in (3.1) is assumed to satisfy: 

w G C^+^n) for some < 7 < 1, gf(-a) = ff (a) for t = 0, 1,2. (3.4) 
We consider a generalized weak solution in the following sense OH]: 
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Definition 3.1. A generalized weak solution of problem (P) is a function h satisfying 



h e C'Jf''^\QT)nL^{0,T;H\Q)), (3.5) 
hteL'iO,T;{H\n)y), (3.6) 

h e CtiiT), ^/Jih) (ao/i... + a^K + a2w') E 1^7), (3.7) 



where 7 = \{h = OUt = 0) and h satisfies (3.1) in the following sense: 

T 

J {ht{-,t),(f)) dt - j j f{h){aoh^^^ + aih^ + a2w'{x))(f)^dxdt 
y 

— hcpxdxdt = (3-8) 



Qt 



for all (j) G C^{Qt) with (p{—a, ■) = 0(a, ■); 



h{-,t) — > /i(-,0) = Hq pointwise & strongly in L^{Q) as t 0, (3.9) 
h{-a, t) = h{a, t) Vt e [0, T] and g(-a, t) = g(a, t) (3.10) 
/or i = 1,3 at all points of the lateral boundary where {h 0}. 



Because the second term of (3.8) has an integral over 7 rather than over Qt, the general- 



ized weak solution is "weaker" than a standard weak solution. Also note that the first term 



of (3.8) uses ht G L^(0,T; (if^(f2))'); this is different from the definition of weak solution 
first introduced by Bernis and Friedman [5]; there, the first term was the integral of kept 
integrated over Qt- 

We first prove the short-time existence of a generalized weak solution and then prove that 
it can have additional regularity. In Section [i] we prove additional control for the norm 
which then allows us to prove long-time existence. 

Theorem 1 (Existence). Let the nonnegative initial data Hq G H^(Q) satisfy 

J^)dx<oo, (3.11) 
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and either 1 ) hQ{—a) = /io(a) = or 2) /io(— a) = hQ{a) ^ and ^|r(— = '^^i'^) holds for i 
1, 2, 3. Then for some time Tioc > there exists a nonnegative generalized weak solution, h, 



on Qtioc the sense of the definition 3.1 . Furthermore, 

heL\0,Tioc;H\n)). 



Let 



^o{T) ■= h {aohl{x, T) — aih'^{x, T) — 2a2w{x)h{x, T)) dx, 



and 



BoiT) := 2 ( ^ + „^ 



(i? I 0-0 U /||2 
II 2 II ^ 



oo, , \\H-,t)\\Uin)dt. 





then the weak solution satisfi,es 



{h>0} 



and 



hlix.Tinr) dx ^ e^«(^'-) / hlJx) dx 



(3.12) 
(3.13) 



^o{Tioc) + jj h\a^h^^^ + ai/i, + a2w'fdxdt ^ £o(0) + KT^o,, (3.14) 



u^y^ju., Hoc) a - j li^^yxjax (3.15) 

n n 
where K = \a2as\ \\w'\\ooC < oo. The time of existence, Tioc, is determined by oq, ai, a2, 
\\w'\\2, \W\\oo, M, f ho, \\ho^\\2, and J I/Hq. 



There is nothing special about the time Tioc in the bounds (3.14) and (3.15); given a 



countable collection of times in [0,T;oc], one can construct a weak solution for which these 
bounds will hold at those times. Also, we note that the analogue of Theorem 4.2 in [S] also 
holds: there exists a nonnegative weak solution with the integral formulation 



{ht{-, t),(j)) dt + ao J J {3h hxhxx(f>x + h hxx4>xx) dxdt 

Qt 



(3.16) 



-Ilia,h. + a2.' + a,h)^^dxdt^0. 



Qi 
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Theorem 2 (Regularity). If the initial data from Theorem^ also satisfies 

J hQ~^{x) dx < oo 
n 

for some — 1/2 < a < 1, a ^ then there exists < T^^J < Tioc such that the nonnegative 

Q a+2 / N 

has the extra regularity /i 2 g L'^{0, TI^^J; H^{Q)) 

andh— eL2(0,T/°j;iy2(r])). 

The solutions from Theorem [2] are often called "strong" solutions in the thin film literature. 
If the initial data satisfies / h'^~^ dx < 00 then the added regularity from Theorem [2] allows 
one to prove the existence of nonnegative solutions with an integral formulation [7] that is 
similar to that of ( |3.16[ ) except that the second integral is replaced by the results of one more 
integration by parts (there are no hxx terms). We also note that if one considered problem 
(P) with nonlinearity f{h) = with < n < 3, then Theorems [T] and [2] would hold for 
general nonnegative initial data Hq e H^{Q); no "finite entropy" assumption would be needed 
[H |3]. Finite entropy conditions (/ /iq~" dx < 00 and / /iq"*"^"" dx < 00) would be needed 
to obtain the results for n > 3. 



3.1 Regularized Problem 

Given 6,e > 0, a regularized parabolic problem, similar to that of Bernis and Friedman [B], 
is considered: 

ht + {f5e{h){aohxxx + aihx + a2w'{x)))^ + a^hx = 0, (3.17) 
-a, t) = £(a, t) for t > 0, ^ = 0, 1, 2, 3, (3.18) 
h{x,0) = ho,5e{x) (3.19) 



d'h . 



where 



fseiz) ■.= fe{z)+5 



\z\+e 



+ 5 WzeR\ 5> 0,6 > 0. 



(3.20) 



The 5 > in (3.20) makes the problem (3.17) regular (i.e. uniformly parabolic). The 
parameter e is an approximating parameter which has the effect of increasing the degeneracy 
from f{h) ~ \h\^ to fe{h) ~ h"^. The nonnegative initial data, Hq, is approximated via 

ho,se = ho^s + £^ e C^^^in) for some < ^ < 2/5 and 7 from ^ 

d'ho. 



i-a) 



(a) for i = 0, 3, 



ho,5£ — > ho strongly in H^{^1) as 5,e 



(3.21) 



0. 
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The e term in (3.21 ) "lifts" the initial data so that it will be positive even if 5 = and the 5 
is involved in smoothing the initial data from H^iVL) to C^^'~^{Vl). 

By Eidelman |T6], Theorem 6.3, p. 302], the regularized problem has a unique classical 
solution hse £ C^y^'^^'^^^{VL x [0, r^e]) for some time r^^ > 0. For any fixed value of 5 and e, 
by Eidelman [TB| Theorem 9.3, p. 316] if one can prove a uniform in time an a priori bound 
\hie{x,t)\ < Ass < oo for some longer time interval [0,T/oc,<5e] (Tioc,Se > Tse) and for all 
X G then Schauder-type interior estimates [IK, Corollary 2, p. 213] imply that the solution 
hse can be continued in time to be in C^'^"''^~^'^^'^{Q x [0,Tioc,Se])- 

Although the solution /i^ is initially positive, there is no guarantee that it will remain 
nonnegative. The goal is to take 5 — > 0, e — > in such a way that 1) Tio^se Tioc > 0, 2) 
the solutions hse converge to a (nonnegative) hmit, h, which is a generahzed weak solution, 
and 3) h inherits certain a priori bounds. This is done by proving various a priori estimates 
for hse that are uniform in 6 and e and hold on a time interval [0, Tioc] that is independent of 
S and e. As a result, {hse} will be a uniformly bounded and equicontinuous (in the Cl^f'^^^ 
norm) family of functions in f2 x [0,T;oc]- Taking 6^0 will result in a family of functions 
{/i^} that are classical, positive, unique solutions to the regularized problem with 5 = 0. 
Taking e ^ will then result in the desired generahzed weak solution h. This last step is 
where the possibility of nonunique weak solutions arise; see |3] for simple examples of how 
such constructions apphed to ht = —{\h\'^hxxx)x can result in two different solutions arising 
from the same initial data. 



3.2 A priori estimates 



Our first task is to derive a priori estimates for classical solutions of (3.17)-(3.21 ) 
lemmas in this section are proved in Section \K\ 

We use an integral quantity based on a function Gss chosen so that 



The 



and Gssiz) > 0. 



(3.22) 



This is analogous to the "entropy" function first introduced by Bernis and Friedman [5]. 
Lemma 3.1. There exists 6q > 0, Sq > 0, and time Tioc > such that if 6 E [0, 5o 



e G (0,5o)? if hse is a classical solution of the problem (3.17)-(3.21) with initial data ho^se, 
and if /lo.fc satisfies (3.21) and is built from a nonnegative function ho that satisfies the 
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hypotheses of Theorem^ then for any T E [0,Tioc] the solution hse satisfies 
I {hlA^,T) + ^ + 25) Gss{hse{x,T))}dx 



(3.23) 



+ ao yy fse{hss)hse,xxxdxdt ^ Ki < 00, 
Qt 



J GseQiseix, T)) dx + oq j j h'j^ ,^^ dxdt < K2 < 00, 



(3.24) 



Qt 



and the energy E,se{t) (see (3.13)) satisfi.es: 



(3.25) 



where = \a2a3\ \\ w'\\ooC < 00. The time Tioc and the constants Ki, K2, Cq, and are 
independent of 5 ande. 

The existence of 5q^ Sq, Tioc, Ki, K2, and K3 is constructive; how to find them and what 
quantities determine them is shown in Section \K\ 

3. l| yields uniform-in-5-and-e bounds for J h'j^^, J Gse{hss), JJ hj^^^, and 



Lemma 



5£,XXX' 



However, these bounds are found in a different manner than in earher 



work for the equation ht = —{\h\'^hx 



for example. Although the inequality (3.24) is 



unchanged, the inequality (3.23) has an extra term involving G^. In the proof, this term 



was introduced to control additional, lower-order terms. This idea of a "blended" ||/ia;||2- 
entropy bound was first introduced by Shishkov and Taranets especially for long-wave stable 
thin film equations with convection [30] . 



Lemma 3.2. Assume Eq and Tioc are from Lemma 



a positive, classical solution of the problem (3.17)-(3.21) with initial data /io,£ satisfying 
Lemma\3.1[ 



3.1. 6 = 0, and e G (0,eo)- U he is 



h'^^(x, T) dx < max 



I 'Il2 



(3.26) 



16 



holds true for all T G [0, Tioc]- Here 

BeiT) ■= 2 



ao 



\he{-MLdt 



The final a priori bound uses the following functions, parametrized by a, 



z :-- 



(3.27) 



(a-l)(a-2) ' {Q-3)(a-2)' V^^^ /.(z)' 

Lemma 3.3. Assume Eq and Ti^c are from Lemma [gj], 5 = 0, and e G (0,£:o)- Assume 



is a positive, classical solution of the problem (3.11)-(3.21) with initial data /io,£ satisfying 
Lemma 3.1. Fix a G (—1/2, 1) with a 7^ 0. If the initial data ho,e is built from Hq which also 
satisfies 

j h'^'^{x) dx<oo (3.28) 



then there exists Sg"'' and T^^^ with < e^q'^ < Eq and < T^^^ < Tioc such that 
{hl,{x,T) + G'f\he{x,T))}dx 

+ [[ [^Kh%, + iK-^hl.] dxdt^K,< 00 



(3.29) 



holds for all T G [0,t/q"^] and some constant K4 that is determined by a, Eq, ao, ai, a2, w' , 
Q and Hq. Here, 

ao if < a < 1, 



«047a - 1/2 < a < 



and 



Furthermore, 



7 



l(l-a) 

a{l—a) 
) 

ao 36 



ao ""''g if < a < 1, 

.il±^^^}(l^ ^/-l/2<a<0. 



a+2 a+2 

he^ eL\0,Tioc;H\n)) and he^ L\Q,Tioc]Wl{n)) 



(3.30) 



with a uniform-in-E bound. 
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The a-entropy, / G''Q'\h) dx, was first introduced for a = — 1/2 in [TO] and an a priori 



bound fike that of Lemma 3^ and regularity results like those of Theorem [2] were found 
simultaneously and independently in |3] and [^. 



3.3 Proof of existence and regularity of solutions 



Bound (3.23) yields uniform L°° control for classical solutions /i^, allowing the time of 
existence Tioc,Se to be taken as Tioc for all 6 G (0,5o) and e G (0,£o)- The existence theory 
starts by constructing a classical solution hse on [0, Tioc] that satisfy the hypotheses of Lemma 
3.1 if 5 G (0,5o) and e G (0,£:o)- The regularizing parameter, 5, is taken to zero and one 



proves that there is a limit and that is a generalized weak solution. One then proves 
additional regularity for h^; specifically that it is strictly positive, classical, and unique. It 
then follows that the a priori bounds given by Lemmas |3.1 and 3^ apply to he- This allows 
us to take the approximating parameter, e, to zero and construct the desired generalized 
weak solution of Theorems [T] and [U 

Lemma 3.4. Assume that the initial data Hq^i. satisfies (g. 21 ) and is built from a nonnegative 
function Hq that satisfies the hypotheses of Theorem [7} Fix 5 = and e G (0, Sq) where Sq 
is from Lemma 3.1. Then there exists a unique, positive, classical solution he on [OjT^oJ of 
problem (Poe), see (3.17)-(3.21), with initial data /loe where Tioc is the time from Lemma 

m 

The proof uses a number of arguments like those presented by Bernis & Friedman [5| and 
we refer to that article as much as possible. 



Proof. Fix e G (0,eo) and assume 6 G (0,5o)- Because Gseiz) > 0, the bound (3.23) yields 
a uniform-in-(5-and-e upper bound on \hse{x,T)\ for {x,T) E Q x [0,Tioc]. As discussed in 
Subsection 3.1, this allows the classical solution hse to be extended from [O^tss] to [0,TioJ- 
By Section 2 of [5], the a priori bound (3.23) on T) ||2 imphes that /i^e G Cl^^t'^^^iQTi^c) 



and that {hse} is a uniformly bounded, equicontinuous family in Qtioc- By the Arzela- 
Ascoli theorem, there is a subsequence {Sk}, so that /i^^.^ converges uniformly to a limit 

he e C'jf''/\Q^J. 

We now argue that is a generahzed weak solution, using methods similar to those of 
[HI Theorem 3.1]. 



By construction, is in Cl^f'^^^iQxioJ^ satisfying the first part of (3.5). The strong 



convergence he{-,t) — > he{-,0) in L'^{Q) follows immediately. The uniform convergence of 



hs^e to he implies the pointwise convergence h{-,t) — * h{-,0) = /iq, and so hg satisfies (3.9) 
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Because hse is a classical solution, 



hs£(j) dxdt 



f5e{h5e){O'0hse,xxx + aihse^x + a2w'{x))(f)x dxdt 



'^3 J J hse4>x dxdt = 0. 

Qt, 

^ loc 



The bound (3.23) yields a uniform bound on 



^ / / dxdt 



(3.31) 



for 5 e (0, So). It follows that 



4 / / {a-ohscxxx + aihse,x + a2w'{x))(f)x dxdt ^0 as 4 0. 



Introducing the notation 

Hse ■= fse{h5e){O'0h5e,xxx + 0'lh5e,x + a2w' (x)) + a^hse 



the integral formulation (3.31) can be written as 



hs£(j)dxdt 



Hs£{x,t)(l)x{x,t) dxdt. 



(3.32) 



(3.33) 



By the L°° control of hs^ and the energy bound (3.25), Hse is uniformly bounded in L'^{Qti^J. 
Taking a further subsequence of {4} yields Hs^e converging weakly to a function if^ in 
L?{Qti„^)- The regularity theory for uniformly parabohc equations implies that /i^e,*, /i5e,x, 
hsE,xx, hse^xxx, and hse,xxxx converge uniformly to h^^t, • • • , hs,xxxx on any compact subset of 



{hs > 0}, implying (3.10) and the first part of (3.7). Note that because the initial data ho,£ 
is in the regularity extends all the way to t = which is excluded in the definition of 7 
in iKTh. 
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The energy ^ssiTioc) is not necessarily positive, 
combined with the L°° control on hse, ensures that £ 



However, the a priori bound (3.23), 
SeiTioc) has a uniform lower bound. As 



a result, the bound (3.25) yields a uniform bound on 



jj fss {h5e)iaoh5e,xxx + aihse,x + a2w'{x) f dxdt. 
Using this, one can argue that for any < a 



(SfjEVX I 



dxdt < Ccy^l'^ 



for some C independent of 5, e, and a. Taking 
conclude 



and using that o is arbitrary, we 



fe (/ie)(ao/ie,xxx + O-lhs^x + a2w' (x)) X, + «3^£- 

{he > 0} 



As a result, taking 5^ in (3.33) imphes /i^ satisfies (3.8) 



The bound (3.23) yields a uniform bound on fj fse{hse)^\ 



Is XXX which can be used in a 



similar manner as above to argue that the second part of (3.7) holds. The bound (3.23) 



also yields a uniform bound on J h'g^^{x,T) dx for every T G [O,r/oc]- As a result, {/i^^e} is 
uniformly bounded in 

L°°(0,r,„,;ifi(fi)). 

Therefore, another refinement of the sequence {5k} yields {hs^e} weakly convergent in this 



space. As a result, G L°°(0, T/oc; -f^^(f^)) and the second part of (3.5) holds. 



Having proven then he is a generalized weak solution, we now prove that is a strictly 
positive, classical, unique solution. This uses the entropy / Gse{hse) and the a priori bound 
(3.24). This bound is, up to the coefficient oq, ident ical t o the a priori bound (4.17) in [5]. 
By construction, the initial data ho,£ is positive (see (3.21)), hence / Ge(/io,e) dx < oo. Also, 
by construction f£{z) ~ z"^. for z <C 1 This implies that the generalized weak solution he is 
strictly positive [5], Theorem 4.1]. Because the initial data /loe is in C'^{il)^ it follows that 

-(4,1 ' 



is a classical solution in C^'t(QT,oJ- This imphes that h^^-.t) — > 
The proof of Theorem 4.1 in [5j then implies that is unique. 



0) strongl}|2 in C\n). 

□ 



^ Unlike the definition of weak solution given in |5], Definition 
converges to the initial data strongly in //^(fi). 



3.1 



does not include that the solution 
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Proof of Theorem\^ As in the proof of Lemma 3.4, following [5], there is a subsequence 



{sk] such that /ig^ converges uniformly to a function h G C].^f'^^^ which is a generahzed weak 
solution in the sense of Definition 3.1 with f{h) = \h\^. 



The initial data is assumed to have finite entropy: / l//io < oo. This, combined with 
f{h) = \h\^, implies that the generalized weak solution h is nonnegative and the set of points 
{h = 0} in Qti^^ has zero measure [SI Theorem 4.1]. 

To prove (|3.14|), start by taking T = T/oc in the a priori bound (|3.25l). As Sk 0, the 



right-hand side of (3.25) is unchanged. First, consider the £fc — > limit of 



£efc(71oc) = \ j ((^ohl^,^{x,Tioc) - aihl^{x,Tioc) - 2a2w{x)he^{x,Tioc)) dx. 



By the uniform convergence of /ig^. to /i, the second and third terms in the energy converge 



strongly as 0. The bound (3.25) yields a uniform bound on { T;oc) dx}. 



Taking a further refinement of {sk}, yields hef.^x{-,Tioc) converging weakly in L^(f2). In a 
Hilbert space, the norm of the weak limit is less than or equal to the liminf of the norms 



< liminf 



loc ) 



of the functions in the sequence, hence £o(^/oc 

on // fs{hs) (ao/i£,xxx + • • • )^ dx also follows from (3.25). Hence \/hJJhk) {aoh 



A uniform bound 
+ ...) 



converges weakly in L'^{Qj'i^J, after taking a further subsequence. It suffices to determine 
the weak limit up to a set of measure zero. Because h > and {h = 0} has measure zero. 



it suffices to determine the weak hmit on {h > 0}. As in the proof of Lemma 3.4, the 
regularity theory for uniformly parabolic equations allows one to argue that the weak limit 
is /i^/^ (aohxxx + ■ ■ ■) on {h > 0}. Using that 1) the norm of the weak limit is less than or 
equal to the lim inf of the norms of the functions in the sequence and that 2) the lim inf of a 



sum is greater than or equal to the sum of the liminfs, results in the desired bound (3.14). 
It follows from (|3.24|) that h^^^^xx converges weakly to some v in L'^{Qti^J, combining with 



strong convergence in L^(0, T; of hg^. to h by Lemma B.l and with the definition 
of weak derivative, we obtain that v = h^x and h E L^{0,Tioc;H^in)) that implies ( |3T2| . 
Hence h^t ^ ht weakly in L^(0,T; (if^(fi))') that implies (3.6). By Lemma B.2 we also have 



3.3 



□ 

The initial data is assumed to have finite entropy 
holds for the approximate solutions {he^} where 



Proof of Theorem^ Fix a G (—1/2, 1' 
jG^^\ho{x)) dx < cx), hence Lemma 
this sequence of approximate solutions is assumed to be the one at the end of the proof of 
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Theorem |T] By ( [33o| 



a+2 



he,^^ > is uniformly bounded in £fc in L^(0, T/oc; 



and 



a+2 



is uniformly bounded in Sk in LF'{{i,Tioc]Wl{Vi)). 



Taking a further subsequence in {e^}, it follows from the proof of fHl Lemma 2.5, p. 330], 



a+2 



these sequences converge weakly in L^(0, T/oc; -f^^(f^)) and L^(0, T/oc; Wi(f^)), to h and 

a+2 

h 4 respectively. □ 



4 Long-time existence of solutions 

Lemma 4.1. Let h G if ^(^2) be a nonnegative function such that J h{x) dx = M > 0. Then 



2 4 



(4.i; 



Note that by taking to be a constant function, one finds that the constant M^/|fi| in 



(4.1 ) is sharp. 



Proof. Letv = h- M/|fi|. By (A.3) 



Hl^n) ^ (DM / v'idx 



If I dx 



Hence 



1 

2 / r \ 3 

/i^ da; 



dx] 



n 



1 

'3^|/ / /,2 \ ^ I M'^ 



hi dx (2M) 3 + 



\n\ ■ 



a 



□ 
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Lemma 4.1 and the bound (3.14) are used to prove control of the generalized weak 



solution constructed in Theorem [H 
Lemma 4.2. Let h be the generalized solution of Theorem\I\ Then 



f\\H;Ti 



ioc)\\m(n) 



<EoiO) + KTi,c + K, 



(4.2) 



where £o(0) is defined in (3.13), M = f Hq, K = \a2a3 



\w 



-,C and 



|a2|||w||oo^ 
|a2||k||ooM + M' 



2|n| 



if aQ + ai < 0, 
otherwise. 



Note that if the evolution is missing either linear or nonlinear advection (02 = or = 
or as = 0) then Lemma 4.2 provides a uniform-in-time upper bound for ||/i(-,T;oc)||hi- 



For the equation (1.3) which models the flow of a thin fllm of hquid on the outside of a 

— |, as = 1, w{x) = sinx, and \Q\ = 2n. In this 



rotating cylinder one has ao = cti 



case, the bound (4.2) becomes 



f , 0-2 



||M.,T,,,)||^i(^) ^ £o(0) + ICTioc + lM + M' (|^+ I) 



where 2£o(0) = /(x/3 {hl^ - hi) + 2/i/3 sin(x) ho )dx. The bound (4.2) actually holds 
true for all times for which h is strictly positive. Recalhng the deflnition (11. ip of x, one sees 



that the control is lost as x ^ (i.e. as a/{upBjjj) — > 0); for example in the zero surface 
tension limit. 



Proof By (3.13) 



j hl{x, T) dx = £o(T) + fj T) dx + a2 J h{x, T) w{x) dx. 
a, n n 



The hnear-in-time bound (3.14) on £o(Tioc) then imphes 



f\\h{;Tioc)\\l.<8.o{0) + KTioc + ^ J h^ dx + \a2\\\w 



00 



M. 



(4.3) 



with K = \a2a2,\\\w'\\ooC. 
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Case 1: + cii < The third term on the right-hand side of (4.3) is nonpositive and can 



be removed. The desired bound (4.2) follows immediately. 



Case 2: ap + ai > By Lemma 4.1 and Young's inequality (A. 5) 



^ I dx< ^ I eiiVft ( f hldx\\^ 



n 



< ^ I hl{x,Tioc) dx + M 
n 



n 

2 f 2^(ao+ai)3/2 qQ+qi 

2\Q\ 



(4.4) 



Using this in (4.3), the desired bound (4.2) follows immediately. 



□ 



This control in time of the generalized solution is now used to extend the short-time 
existence result of Theorem [T] to a long-time existence result: 

Theorem 3. Let Tg be an arbitrary positive finite number. The generalized weak solution h 
of Theorem^ can be continued in time from [0,Tioc\ to [0,Tg] in such a way that h is also 
a generalized weak solution and satisfies all the bounds of Theorem^ (with Tioc replaced by 

Similarly, the short-time existence of strong solutions (see Theorem [2]) can be extended 
to a long-time existence. 

Proof. To construct a weak solution up to time Tg, one applies the local existence theory 
iteratively, taking the solution at the final time of the current time interval as initial data 
for the next time interval. 
Introduce the times 



N-l 



= To < Ti < T2 < ■■■ < Tiv < ... where T^:=^T„ 



loc 



(4.5) 



n=0 



and Tn^ioc is the interval of existence (A. 20) for a solution with initial data h{-,Tn): 
^n,Zoc:=iLniin<[l, | / hl{x,Tn) + 2f^Go{h{x,Tn)) dx 



(4.6) 
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The proof proceeds by contradiction. Assume there exists initial data Hq, satisfying the 
hypotheses of Theorem [TJ that results in a weak solution that cannot be extended arbitrarily 
in time: 



ETn,loc = T* < OO =^ lim Tn^ioc = 0. 
n— >oo 

k=0 

From the definition ( |4.6[ ) of Tn^ioc, this imphes 

lim [ {hl{x,Tn) + 2fGo{h{x,Tn))) dx = oo. 



(4.7) 



By (4.2), 



By (3.14), 



Combining these. 



4 



^ j hl{x,T^) dx < £o(T„_i) + KT^_,^i,, + 
n 

So(^n-l) ^ £o(^n-2) + K T„_2,/oc- 



hl{x, Tn) dx < £o(^n-2) + K {Tn-2,loc + Tn-l,loc) + ^3- 



Continuing in this way. 



J J hl{x,T^) dx < £o(0) + KT^ + K, (4.8) 
n 

By assumption, T„ — ^ T* < oo as n ^ oo hence / h'^{x, T„) remains bounded. Assumption 



(4.7) then implies that / Go{h{x,Tn)) — oo as n — ^ oo. 



analogue of T„ and Tn^ioc,e, defined by (lA.lQh, be the analogue of T„^/oc By (|A.16l), 



To continue the argument, we step back to the approximate solutions /i^. Let T„£ be the 

(4.9) 



J Ge{he{x,Tn,e)) dx < j Ge{he{x,Tn-l,e)) dx 



+ Cg j max < 1, y /i^ ,^(x, T) (ix > dT 
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Using the bound (3.25), one can prove the analogue of Lemma 4.2 for the approximate 



solution he- However the bound (4.2) would be replaced by a bound on T)||/i^i which 



holds for all T G [0, Teniae]- This bound would then be used to prove a bound like (4.8) to 
prove hnear-in-time control of / /i^^(x,T) for all T G [0,T„ e]. Using this bound, 



hl^{x, T) dxdT < 



ao 



{ZM + KT + K^) dT 



4_ 

ao 



[£,(0) +^3+1 {Tn-l,e + Tn,e)] T^-l,loc,e- 



(4.10) 



If the initial data is such that 4/ao (£e(0) + ^3) < 1 then before using (4.10) in (4.9) we 



replace by a larger value so that 4/ao (£e(0) + K^) > 1. Using (4.10) in (4.9), it follows 
that 



Ge{hs{x,Tn,e)) dx 



[4.11] 



< j Geihei XyTn—l^e)) dx + (ci + f3(Tn—l,e ~l" T^^e)) Tn—l^loc,e 



for some a and /3 which are fixed values that depend on the coefficients of the PDE, and 
(possibly) on the initial data hQ ^- 

One now takes the sequence {sk} that was used to construct the weak solution of Theorem 
[l|on the interval [T„_i,T„]. Taking ek 0, ( |4.1lD yields 



Go{h{x, T„)) dx < J Goih{x, T„_i)) dx + (a + /5(T„„i + T„)) Tn-i^ioc- 
n a 



(4.12) 
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Applying (14. 121) iteratively. 



J Go{h{x,Tn)) dx < J Go{h{x))dx + J2Tk,ioc {a + [3 {Tk + Tk+i)) 

/n—l 
Go{ho{x)) dx + J2Tk,ioc {a + 2(3T*) 
1 n 



Go{ho{x)) dx+{a + 2pT*)Tn. 



This upper bound proves that / Go{h{x,Tn)) cannot diverge to infinity as n — > oo, finishing 
the proof. 

□ 



Under certain conditions, a bound closely related to (4.2) implies that if the solution of 



Theorem [T] is initially constant then it will remain constant for all time: 



Theorem 4. Assume the coefficients ai and 02 in {1.8) satisfy ai > 0, 02 = and \Vt\ < 
4ao/|a.i|. If the initial data is constant, hQ = G > 0, then the solution of Theorem^ satisfies 
h{x, t) = G for all X E Cl and all t > 0. 



The hypotheses of Theorem [4] correspond to: the equation is long-wave unstable (ai > 0), 
there is no nonlinear advection (02 = 0), and the domain is not "too large". 

Proof. Consider the approximate solution h^. The definition of £-e{T) combined with the 



linear-in-time bound (3.25) implies 



f j h%{x,T) dx < £,(0) + irT+ ^ j hi dx+\a2\ 



\w\\ooMe 



(4.13) 



where = J ho^^dx. Applying Poincare's inequality (A. 2) to = h^ — Me/\VL\ and using 
jhldx = jvldx + M2/|fi| yields 



2 8 



j hl,{x,t) dx < 6,(0) + KT.joc + ^ + \a2\ 



\w\\ooMe. 
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If /io,e = = C + and a2 = (hence K = 0) this becomes 

(f - ^) J hUx,T) dx < {a, - \a,\)^. 
n 

If ai > and \Q\ < Aao/ai then / h'^ ,^{x,T) dx = for all T G [0,Te^/oc] and that this, 
combined with the continuity in space and time of h^, implies that /i^ = on Qt^ 

Taking the sequence {sk} that yields convergence to the solution h of Theorem [l| h = C 
on Qt,„,. 

□ 



5 Strong positivity of solutions 



Proposition 5.1. Assume the initial data ho satisfies ho{x) > for all x & u C Q where u 
is an open interval. Then the weak solution h from Theorem [I] satisfies: 

1) h{x,T) > for almost every x E u, for all T G [0,Tioc] 

2) h{x, T) > for all x E u, for almost every T G [0, T^oc] 



The proof of Proposition 5.1 depends on a local version of the a priori bound (3.24) of 
Lemma |3. II 



Lemma 5.1. Let u ^VL be an open interval and ( G C^(f2) such that ( > onu, supp ( - 
and (C^)' = on dQ. If u = Q, choose ( such that C{~'^) = C('^) > 0. Let ^ := 

If the initial data Hq and the time Tioc are as in Theorem^ then for all T G [OjT^oJ 
weak solution h from Theorem [I] satisfies 



the 



h{x,T) 



dx < oo 



(5.i; 



The proof of Lemma [5T] is given in Appendix|X] The proof of Proposition 5T is essentially 
a combination of the proofs of Corollary 4.5 and Theorem 6.1 in [S] and is provided here for 
the reader's convenience. 



Proof of Proposition 
Lemma 5.1 Hence, 



5.1 



5.1) 



Choose the localizing function ({x) to satisfy the hypotheses of 
holds for every T G [0,Tioc]. 
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First, we prove h{x, T) > for almost every x G u;, for all T G [0, Tioc]- Assume not. Then 
there is a time T G [0, TjoJ such that the set {x \ h{x^ T) = 0} fl u; has positive measure. 
Then 

oo> j ax)j^^dx> J ^ix)j^^dx = oo. 

n {h{-,T)=o}nuj 

This contradiction imphes there can be no time at which h vanishes on a set of positive 
measure in u, as desired. 

Now, we prove h{x,T) > for all x G u;, for almost every T G [OjT^oc]- By (3.12), 
h^^{-,T) G L^{Vt) for almost all T G [0,Tio,] hence h{-,T) G C^/^{Vt) for almost all T G 
[0,Tio,]. Assume Tq is such that /i(-,To) G C^/'^{Vl) and /i(xo,ro) = at some Xq G ciJ. Then 
there is a L such that 

h{x,To) = \h{x,To) - h{xo,To)\ < L\x-Xo\^^'^. 

Hence 

^(^)m^ rfx > W C(x)|x - xol"^^^ dx 



This contradiction imphes there can be no point xq such that h{xo,To) = 0, as desired. Note 
that we used ,^ > on a; and Xq G to conclude that the integral diverges. 

□ 



A Proofs of A Priori Estimates 

The first observation is that the periodic boundary conditions imply that classical solutions 



of equation (3.17) conserve mass: 



J hs£{x,t) dx = J ho^se{x) = < oo for all t > 0. (A.l) 



Further, (3.21) implies Mse ^ M = J /iq as e, 5 — > 0. The initial data in this article have 
M > 0, hence Mss > for 5 and e sufficiently small. 

Also, we will relate the norm of h to the norm of its zero-mean part as follows: 



|Mx)|< h{x)-^^ 



M_ 



l<2^-' \\v\\l+[ 
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where v := h — M/\fl\ and we have assumed that M > 0. We will use the Poincare inequality 
which holds for any zero-mean function in H^{Q) 



\vE<bi\\v,\\l l<p<oo 



(A.2) 



with bi = \n\P/{p2P-^). 

Also used will be an interpolation inequality [2T1 Th. 2.2, p. 62] for functions of zero 
mean in H^{Q): 



where r > 1, p > r, 



\vE < 62 \\vx 



l/r-l/p 



I 'JP 



l/r+1/2' ^ 

It follows that for any zero-mean function v in H^{Q) 



(l + r/2) 



ap 



hE<b4h,\\l + b,Ml 



(A.3) 



(A.4) 



where 



^ l^|(2-p)/p if l<p<2 



b^r'^" b. 



if 2 < p < 00 ' 



2P-^ b 



'3, 



P-i 



To see that (A.4) holds, consider two cases. If 1 < p < 2, then by (A.2), ||f ||p is controlled by 
||v3;||p. By the Holder inequality, \\vx\\p is then controlled by ||fx||2- If p > 2 then by (A.3), 
ll^llp is controlled by ||i'a;||2ll'^ll2~" where a = 1/2 — 1/p. By the Poincare inequality, ||f ||2~" 
is controlled by ||va;|l2~"- 

The Cauchy inequality ab ^ ea^ + b"^ /{At) with e > will be used often as will Young's 
inequality 



ab < eaF + 



q {ep)i/P ' 



p q 



1, e > 0. 



(A.5) 



Proof of Lemma \3.1\ In the following, we denote the classical solution hse by h whenever 
there is no chance of confusion. 

To prove the bound (3.23) one starts by multiplying (3.17) by —hxx, integrating over Qt, 
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and using the periodic boundary conditions (3.18) yields 



2 , h^{x,T) dx + ao j j fse{h)h^^^dxdt 
n Qt 



I j hQ^5e,x^{x)dx - ai jj f^{h)h^h^^^dxdt + Sai jj hl^ dxdt 

Q Qt Qt 



T VT 

- a, jj &(A)u/ft„ dMt - Sa, jj ^ h„. d.dt. 

Qt Qt 

The periodic boundary conditions were used above to conclude 



(A.6) 



as / / h^h^^ dxdt = y 1 1 {hl)xdxdt = f I hl{a,t) - hl{-a,t)dt = 0. 



Qt 



Qt 



The Cauchy inequality is used to bound some terms on the right-hand of (A.6): 



ai / / fse{h)KK^^ dxdt ^ f / / fseWhl^^ dxdt + ^ / / fseWhl dxdt, (A.7) 



"2 Jj f5eih)w'Kxxdxdt ^ ^ jj fseih)hl^^dxdt + ^ Jj fseih)w'^ dxdt, 

Qt Qt 



Qj 



^0-2 J J w'hxxx dxdt < 5 

Qt Qt 



f / / Kxx dxdt + S^^T I dx. 



(A.8) 
(A.9) 
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Using ([0)-(|X9l) in ([X6| yields 



\ J hl{x,T)dx + f J J fse{h)hl,,dxdt 
n Qt 



\ j h,se/dx+'^^ J J fs,{h)hldxdt+f^ J J fseih)w"dxdt 
n Qt Qt 

+ 6 ai / / h^^ dxdt + 6 T w'^ dx 



Qt 



J ho,Se/dx+f^ jj Ihl^'hldxdt+^^Jw'Wl, jj \hfdxdt 



Qt 



Qj 



+ (5 oi / / h^^ dxdt + S T I w dx 



Qt 



a 



(A.IO) 



(A.11) 



Above, we used the bound fse{z) < \z\^. By the Cauchy inequahty, bound (A.4), and bound 



(A.3) 



\h\^hl dxdt < I dxdt + k dxdt 



2 11 I 2 

Qt Qt Qt 

^'i j [j hl'^A dt+'-iMlT+f j \\h,A-M2\\h.{;t)\\ldt 
\n / 

~^ II ^^^^ + ci J i J hi dx\ dt + C2T 



(A.12) 



Qt 



\n 



where ci = 6^/8 + h/2 and ca = 65/2. 
By (itg, 

T 



3/2 



\h\^ dxdt I \ I hldx\ dt + b^ Ml T. 



Qt 



\n 



(A.13) 
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From (lA.lip, due to (|A.12p-(|A.13p, we arrive at 



\ j hlix,T) dx + fJJ fseih)hl,,dxdt 
n Qt 



^ I y ^o,(5e,a:^ dx + C3 J J h^.^ dxdt + C4 J i J hi dx\ 
n Qt \n / 



dt 



T / \ 3/2 

+ C5 J i J hi dx \ dt + cqT 
\n 



— i y ^o,5e,x'^ dx + C3 J J hl^ dxdt + '^7 y ^^-^ I ' ( y '^^ 



dt 



where 



0.0 



ao I 



C6 = + ^l|u;'||L&5M|, + d'^Jw'Wl cr = c, + c, + c,. 



(A.14) 



Now, multiplying (3.17) by G'g^{h), integrating over Qt, and using the periodic boundary 
conditions (|3.18l), we obtain 



Gs£{h{x,T)) dx + ao J J h^^dxdt = J Gseiho^se) dx + ai J J h^dxdt 

n Qt ^ Qt 

— as jj {Gse{h))xdxdt + a2 jj w'h^dxdt. (A. 15) 

Qt Qt 

By the periodic boundary conditions, 



iGssih))xdxdt = J {Gseih{a,t))-Gssihi-a,t)))dt = 0. 

Qt 
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Thus, from (fATsj), we deduce 



Gs£{h{x,T)) dx + ao J J h^^dxdt 
n Qt 

^ J Gseiho^ss) dx + ai J J h^dxdt + \a2\\\w'\\2 
n Qt on 

T 




1 

dx 1 dt 



where 



^ J Gseiho^ss) dx + cs J max 1, / h^{x,t) dx ^ dt, 
a, V n 

cs = + |a2|||w'||2- 



Further, from (A. 14) and (A. 16) we find 



h^dx + ^ J Gse{h{x, T)) dx + ao j j f5eih)hl^^ dxdt 



Qt 



dxdt 



< J ho,5e/ dx+^^ij GssiHx, T)) dx + ao jj . 
n \n Qt 

+ 2c7 J max < 1, I y hl{x,t) dx\ \ dt < J hQ^se,x^ dx 
I \n / J n 

J Gseiho^se) dx + Cs J laax l^^J ^xi^^'^) dx > dt 



+ 



2£3 
ao 



K n 

3- 



+ 2c7 J max < 1, | y hl{x,t) dx\ \ dt < J /lo.fe.a:^ dx 
I \n / j n 

+ ^ J Gssiho^se) dx + cg J max < 1, | y h1{x,t) dx\ \ dt 
n K \n / ) 



(A.16) 



(A.17) 
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where Cg = 2c^c^/aQ + 2c7. 

Applying the nonlinear Gronwall lemma [12] to 



J. 

v{T) < f (0) + Cg y" max{l, v^{t)} dt 



with v{t) = J hl{x,t) + 2c3/ao Gseih{x^t)) dx yields 

'v{0) + Cgt if t<to 



From this, 



v{t) < { |^(l-2c9(t-to))~'^' if t>to 



hl{x,t)+2f^Gse{h{x,t))dx 



^9 if v{0) < 1 
if v{0) > 1. 



(A.18) 



< A/2max|l,y" {ho,se/{x) + 2'^Gss{ho,se{x))) dx^ = Kse < oo 
for all t G [0, Tse,ioc\ where 

4 min J 1, I y (/io,fe,x'(a:) + 2f^Gse{ho,se{^))) dx 



S£,loc ■ 4(,g 



(A.19) 



Using the 5 — > 0,5 ^ convergence of the initial data and the choice oi 6 e (0,2/5) (see 



(3.21)) as well as the assumption that the initial data ho has finite entropy (3.11), the times 



Tse,ioc converge to a positive limit and the upper bound K in (A.18) can be taken finite and 



independent of 6 and e for 6 and e sufficiently small. (We refer the reader to the end of the 
proof of Lemma 5.1 in this Appendix for a fuller explanation of a similar case.) Therefore 



there exists 5o > and > and K such that the bound (A.18) holds for all < 5 < (5o 
and < e < Eo with K replacing Kse and for all 



< t < T, 



loc 



_9_ 
10 



lim 



Se.loc- 



(A.20) 
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Using the uniform bound on f h?^ that (A. 18) provides, one can find a uniform-in-5-and-e 



bound for the right-hand-side of (A. 17) yielding the desired a priori bound (3.23). Similarly. 



one can find a uniform-in-5-and-e bound for the right-hand-side of (A. 16) yielding the desired 



a priori bound (3.24). 



To prove the bound (3.25), multiply (3.17) by —a^h^x — o-ih — a2W, integrate over Qt, in- 



tegrate by parts, use the periodic boundary conditions (3.18), and use the mass conservation 



(see (A.l)) to find 



£fe(T) + J J fse{h){aohxxx + aihx + a2w'{x)) dxdt 

Qt 

= £e(0) + 0203 / / hxw{x) dxdt 



= £e(0) — 02^3 J J hw'{x) dxdt 

Qt 

< £2(0) + 10203! II II 00 J J \h{x, t) \ dxdt. 

Qt 



(A.21) 



By the embedding theorem (A. 4) and the bound (3.23), one has 



\h{x,t) \ dxdt < \n\ 



\ n 



\hx{x,t)\^ dx dt + MssT 



< [\n\'^^/K~l + 2MJT, (A.22) 
where Sq and Eq have been taken smaller, if necessary, to ensure that < 2M. Substituting 



(A.22) into (A.21) yields the desired bound (3.25) with the constant 

|a2a3||k'||oo(|^^rv^ + 2M). 



K3 . . , _ 

The parameters 60 and Eq are determined by ao, ai, 02, ||w'||2, ||'U^'||oo, |^|, / ^0, ll^oxlh, 
J l/ho, by how quickly Mse converges to M, and by how quickly the approximate initial 
data (3.21), /?,o,fe, converge to ho in H^{il). 



The time Tioc and the constants Ki, K2, and are determined by 5q, Eq, a^, oi, 02, 
k'lh, Ik'lloo, l^^l, / ll/ioxlh, and / l/Zio- □ 
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Proof of Lemma [X^ In the following, we denote the positive, classical solution by h 
whenever there is no chance of confusion. 



Taking 5 ^ in (AAO) yields 



I hldx + f f,{h)hl,,dxdt 



n 



I I hl^, dx+'l^ll fs{h)hl dxdt +ijjfe 
n Qt Qt 



{h)w' (x) dxdt 



(A.23) 



< i J hl^,dx + J \\hi.,t)\\l J ^^hl{x,t) + dx I dt 
n \ Q 



J hl^,dx+^^ j\\hi;t)\\l max|||^'||2,| hl{x,t) dx\ dt (A.24) 
n K n ) 



Applying the nonlinear Gronwall lemma |T2] to 



v{T) < v{0) + / B{t) max{A,v{t) dt 



with v{t) = J hl{x,t) dx, B{t) = J (2(a? + a^)/ao \\h{- 



Ml] 



dx and A = \\w'\\2 yields 



v{T) < max{A,v{0)} e" 



J Bit) dt 



□ 



Proof of Lemma \3^ In the following, we denote the positive, classical solution hg by h 
whenever there is no chance of confusion. 

Multiplying (3.17) by {G'f\h)y, integrating over Qt, taking 5 — > 0, and using the periodic 
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boundary conditions (iS-lSh, yields 



j G'f\h{x,T))dx + ao jj h^hl^ dxdt + ao^^^ jj h^-^ht dxdt 



(A.25) 



Qt 



j V) dx + aijj h'^hl dxdt Jj h^^^w" dxdt. 

n Qt Qt 



Case 1: < a < 1. The coefficient multiplying ///i" ^/i^ in (A.25) is po sitive and can 



therefore be used to control the term Jf h°'hl on the right-hand side of (A.25). Specifically, 
using the Cauchy-Schwartz inequality and the Cauchy inequality, 



ai y / h°'h^ dxdt ^ 

Qt 



aoa(l— a) 



h"-^ht dxdt + 



3ai 



2aoa{l—a) 



dxdt. 



Qt 



Using the bound (A. 26) in (A.25) yields 



G^^\h{x,T))dx + ao jj h'^hl^dxdt + ao^^^ jj h'^'^ht dxdt 

Qt Qt 

< j G^^\h,,) dx + j I ^^^^ + \^,\\w"U j I h''^' dxdt. 



Qi 



By (A.4) 



2^2 



dxdt ^ b. 



4 I 1 I ■-x 

\n 



hi dx\ dt + hM^+^ T, 



r+i 



^^^t I \ I hldx] dt + hX^^ T. 



Qt 



\n 



(A.26) 



(A.27) 
(A.28) 



(A.29) 



(A.30) 
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Using ([X29| and (|A3o| in (|X27| yields 



Gi"\h{x,T))dx + ao // h"hl^dxdt 



Qt 



+ ao 



q(1— a) 



Jj h'^-^hldxdt< j G^^\ho,)dx 

Qt n 
T / \ T / 

+ di j \j hldx\ dt + d2 J IJ hldx 
\n I \n 

< y G["^(/ioe) dx + di J max < 1, | y /i^ 



2^2 



r+1 



dt + d^T 



dt 



(A.31) 



where 



cii = 64 



3af 



2aoa(l— a) ' 



1+a 



^3 = &5 (2^^ M^^ + »^ M^i) , d, = d,+d, + d. 



Taking 5 in (A. 11) yields 



hi dx + ao JJ fs{h)hl^^dxdt 

Qt 



i^hldxdt + ^\\w'\\l^ jj h^dxdt 

n Qt Qt 

Applying the Cauchy-Schwartz inequality and the Cauchy inequality, 



(A.32) 



2al 
ao 



h^hl dxdt 



Qt 



< 



aoa(l—a) 



h'^-^ht dxdt + 



aga{l-a) 



/i^-" dxdt 



Qt 



Qt 



(A.33) 
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By dXl, 



h^-" dxdt ^ bA 



hi dx\ dt + hM^-" T, 



\n 

T 



dxdt \ hldx \ dt + h^M^ T. 



\n 



Using ( [X33| , ( |X34| and ( [X35| ) in ( [X32| ) yields 

hldx + ao JJ f!r{h)hl^^dxdt ^ J hl^^^dx 



Qt 



+ 



6 



aoa(l-a) / / r a-2T4 



h" dxdt + (is 



dt 



\n 



+ da 



h'^ dx ] dt + dj T < I hl^ ^ dx 



\n 



where 



Add 



+ 



6 



aoa(l-a) / / 7 0-2^4 



/i" dxdt + ds max < 1, dx 



dt 



Qt 



aga(l— a) 



ao 



d? = h 



loo e / ' 



(is = rfs + (^6 + 



G'f\h{x,T))dx 



to both sides of (A. 36) and add 



ao h"h? dxdt 



Qt 
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to the right-hand side of the resulting inequahty. Using (|A.31l) yields 



n 



hl{x, T)dx + j G^^\h{x, T)) dx + aojj fe{h)hl^^ dxdt 
n Qt 

— J ^os,xd'X + J G^^\hQf.) dx + di J max < 1, | y h"^ dx 



(A.37) 



+ c?8 y max < 1, I / hi dx 
I \n 



dt 



< J hl^,^dx + J G^f\hoe) dx + dg j max <j 1, I J hi dx 



where dg = d^ + ds. 

Applying the nonlinear Gronwall lemma [I2] to 



1 

v{T) < v{Q) + dg j max{l, v^-'^''^{t)} dt 



with v{T) = Jhl{x,T) + G^^\h{x,T)) dx yields 

'v{<d) + dgt if T<To:=i^ 



viT) < < 



{l-^dg{T-To)) if T>Tn 



6— a 



V{0)- 2 - ^dgT 



2 

6— « 



if v{0) < 1 
if v{0) > 1 
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From this, 



for all 



{hl{x,T) + G^f\h{x,T))) dx 



(A.38) 



< 46-a max 1 1, y iho/^{x) + G'(")( V(x))) dx^ = K, < oo 



6— g 
2 



< T < T 



n 



mm{Tioc, T^ilc} where Tioc is from Lemma 



The bound (A.38) holds for all < £ < Eq where Eq is from Lemma 3.1 and for all t < 



3.1 



Using the 5 — 0,e — > convergence of the initial data and the choice of 6* G (0,2/5) 



(see (3.21)) as well as the assumption that the initial data has fini te a-e ntropy (3.28) 



the times T^°^l^ converge to a positive limit and the upper bound in (|A.38|) can be taken 

in 



Lemma 5.1 



finite and independent of e. (We refer the reader to the end of the proof of 
this Appendix for a fuller explanation of a similar case.) Therefore there exists e^^^ and K 



such that the bound (|A.38[) holds for all < e < e^^^ with K replacing and for all 



< t < T, 



(a) 



loc 



min < T; 



(A.39) 



where Tj^r is the time from Lemma 



3.1 



Also, without loss of generality, e'"^^ can be taken 
jemma 



3.1 



to be less than or equal to the Eq from 

Using the uniform bound on / h"^ that (A.38) provides, one can find a uniform-in-e bound 
for the right-hand-side of ( A.31[ ) yielding the desired bound 



j G^^\hix, T)) dx + aoJJ h^hl^ dxdt + ao^^^ j j h^'^hl dxdt < Ki (A.40) 



which holds for all < e < ej,"^ and all < T < t/„"\ 
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It remains to argue that (A.40) implies that for all < e < ^o"-* that h"^"^^^ and h'^^'^^'^^'^ 
are contained in balls in L^(0, T; H'^{VL)) and L^(0, T; Wl{Vt)) respectively. It suffices to show 
that 

j j {K'^^^)l^ dxdt <K, jj dxdt < K 



for some K that is independent of e and T. The integral // ^ hnear combination 
of // h"~'^hl, JJ h°'~^hlhxxy and // h'^hl^. Integration by parts and the periodic boundary 
conditions imply 

^ h'^-^ht dxdt = !l K'-^hlK^ dxdt [kAl] 



l-g 
3 



Hence JJ{h. 



a/2+lN2 



is a hnear combination of JJ h°' ^/i^, and JJ h"hl^. By (A.40), the 



two integrals are uniformly bounded independent of e and T hence JJ{h£^'^^^)l^ is as weU, 



yielding the first part of (3.30). 

The uniform bound of JJ{he^^^^^'^)'^ follows immediately from the uniform bound of 



JJ h" "^h^, yielding the second part of (3.30). 



Case 2: — | < a < 0. For a < the coefficient multiplying ^/^^ in (A. 25) is 



negative. However, we wiU show that if a > —1/2 then one can replace this coefficient with 
a positive coefficient while also controlling the term JJ h'^hl on the right-hand side of (A.25). 
Applying the Cauchy-Schwartz inequality to the right-hand side of (A.41), dividing by 



// h°'~'^h'^, and squaring both sides of the resulting inequality yields 



h'^-^h't dxdt < 



h^hl^ dxdt 



Qt 



Qt 



Va < 1. 



(A.42) 



Using (A.42) in (A.25) yields 



Gi''\h{x,T))dx + ao^_ 



l+2a / / T,Q/,2 



h^h^^ dxdt 



(A.43) 



< j G^^\hoe) dx + ai jj h^hl dxdt + ^||w"||oo j I dxdt. 

^ Qt Qt 
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Note that if a > —1/2 then all the terms on the left-hand side of (A. 43) are positive. We 
now control the term J J h^'hl on the right-hand side of (A.43). 
By integration by parts and the periodic boundary conditions 



h^'hl dxdt = -j^ 1 1 h'^+^K^, dxdt 



(A.44) 



Applying the Cauchy-Schwartz inequality and the Cauchy inequality to (A.44) yields 



a,jj h-hl dxdt < ^|[i±^ j j h-hl dxdt + 1 1 h"^' dxdt (A.45) 



Using inequality (A.45) in (A.43) yields 



Gi''\h{x,T)) dx + ao^ff^^ jj h"hl^ dxdt 
n Qt 

< j G^^\hoe)dx + 



(A.46) 



2ao(l+2a)(l+a)2 / / axai ^ ^^^\\UJ ||oo 



dxdt. 



Adding 



ao(l+2a)(l-a) / / 4 
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h'^-'K dxdt 



to both sides of (A.46) and using the inequality (A. 42) yields 



G(")(/i(x,T))cix + ao^_„) 

n Qt 
ao(l+2Q)(l-a) / / La-27 4 



K^h^^ dxdt 



+ 
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j I h'^-^hl dxdt < j G^^\hoe) dx 



(A.47) 



+ 



2.0(112^71+.)- jj h-^'dxdt+j^Jw"\\^ llh^^Uxdt. 

Qt Qi 
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Using ([X29| and (|A3o| in (|X47| yields 



Gi"\hix,T))dx+^-^^ jjh-hl^dxdt 

Qt 

jj K'-^hldxdtK j Gf\ho,)dx 



aa{l+2a)(l-a) 
36 



Qt ^ 

a 



+ ei 



1 nr 
hi dx \ dt + 62 



T / \ 2 

2 



2^2 



hi dx I dt + esT 



\n 



\n 



64 J max < 1, I y hi dx 



dt 



(A.48) 



where 



1*^21 ||„,,"1 



^1 ~ 2ao(l+2a)(l+a)2 "4, 62 - l|oo O4, 



and 64 = ei + 62 + 63. 



Recall the bound (A.32): 



hldx + ao jj fe{h)hl^^ dxdt 



Qt 

< I hl^^dx+'-§ [[h'hldxdt+'-^\\w'\\l II hUxdt. 



n Qt Qt 

As before, by the Cauchy-Schwartz inequality and the Cauchy inequality, 

'-^jj h'hl dxdt < ""(^+^3^^-°) j I' h^'^hl dxdt 
Qt Qt 



+ 



36af 



a3(l+2a)(l-a) 



(A.49) 



(A.50) 
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Using (IXSOl, (IXsl, and ([X35| in (1X491 yields 



Qt n 



+ 



ao(l+2a)(l-a) 
36 



/i^ dx 



dt 



Qt 



\n 



+ e6 j i J hldx\ dt + e7T< j hl^^^ dx 



+ 




ao(l+2a)(l-a) 



(A.51) 



where 



Add 



65 



36af 7 
1 0. 



a^(l+2a)(l-a) 



4, 



2a| ||^,,/||2 

ao 



7 / 36af 

67 = 6,- ' 



'5 I ag(l+2a)(l-.) + H^'H^ ) ' ^8 = 65 + 65 + 67. 



G(")(/i(x,T))da; 



to both sides of (A.51) and add 



ao(l+2a) 
4(l-a) 



h'^hl^ dxdt 



Qt 



to the right-hand side of the resulting inequality. Just as (A. 31) and (A.32) yielded (A.37) 
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(A. 48) combined with the above inequality yields 



hl{x,T)dx+ / Gi''\h{x,T))dx + ao // f,{h)hl^^dxdt 



Qt 



— J ^0£,xd'X + J G'f\hQe) dx + 64, J max < i J ^1 dx 



no. 

{A 
(l \ ' 

lA I h^. dx \ \ dt 

T 



eg j max < 1, | y h^. dx 



(A.52) 



where eg = e4 + eg. 

The rest of the proof now continues as in the < a < 1 case. Specifically, one finds a 
bound 



{hl{x,T) + G^^\h{x,T))) dx 



(A.53) 



< 46-« max <( 1, / ( (a; 



j ( + G^^\ho,e{x))) dx\=K,<CX> 



for all 



6~a 
2 



< T < T 



(") ._ 1 



£,l0C • 69(6-0) 



mm < 1, I j iho,e/ix) + Gi''\ho,,{x))) dx 



The time T^^J is defined as in (|A.39[) and the uniform bound (|A.53[) used to bound the right 
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hand side of (|A.48l) yields the desired bound 



Qt 



+ "0(l+2a)(l-a) / / ^.-2^4 ^^^^ < 



(A.54) 



□ 



Proof of Lemma \5A In the following, we denote the positive, classical solution con- 
structed in Lemma 3.4 by h (whenever there is no chance of confusion). 



Recall the entropy function Gsei^z) defined by (3.22). Multiplying (3.17) by ^{x)G'^^{hse)-, 
taking 5 — > 0, and integrating over Qt yields 

^{x)Ge{h{x,T)) dx — I ^{x)Ge{hQ^i;) dx = —as H ^{x)G'^{h)hx dxdt 

U Qt 

+ II feih){aoh^^x + aih^ + a2w'){^'G'^{h) + ^G'^{h)hx) dxdt 

Qt 

as II i'Ge{h)dxdt + I I ^' fs{h)G'^{h){aoh^x^ + aih^ + a2w') dxdt 

Qt Qt 

+ ^K{aoKxx + aihx + a2w') dxdt =: h + h + h- (A.55) 
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We now bound the terms lo and h. First, 



h --^oJJ e fe{h)G'^{h)K, dxdt -aoJJ r (1 + ah)G'^{h))KK, dxdt 

Qt Qt 

-ai JJ ^"Fe{h) dxdt + 0.2 jj fe{h)G'^{h)w'{x) dxdt 



Qt Qi 



--^oJJ efe{h)G',{h)h,,dxdt + f J J eiildxdt 

Qt Qt 

-ooJJ a'e{.h)G',{h)h,h,^dxdt -oiJJ eFe{h) dxdt 

Qt Qt 

+ 02 J J efe{h)G'^{h)w'{x) dxdt, 

Qt 

where F,{z) := j /,(s)G^(s) ds = -z^A + e/6z- \n{z + e). 


One easily finds that for all £ > and all 2; > 

\Mz)G',(z)\ = ^ Iz, \mG',(z)\ = i (^-+g(y3.) ^ 2. 

To bound a limit on the possible values of e is assumed. Specifically, if < £ < 

(^33 - 3)/4 then for all ^ > 

\FM\<b" + l 

Using these bounds, and recalhng ^ = C^, we bound \I2\: 



\I2\ ^ f 



JJ [l2CC + K''\Cxx\\h\K^\dxdt + 8ao JJ C^\Cx\\ha^\\h^:c\ dxdt 

Qt Qt 

+ ^aoJJ [3C'C' + C'|Cx|] hldxdt+\a,\ J J [6C'C' + 2C'|Ca:x|] dxdt 



Qt Qt 



+ 2\a2\\\w'\\oo JJ C^\Cx\h dxdt + lla^l JJ \i"\dxdt (A.56) 
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for all < e < (v33 — 3)/4. The first two integrals on the right hand side of (A. 56) are 
bounded via the Cauchy-Schwartz inequality followed by the Cauchy inequality yielding: 



(A.57) 
(A.58) 
(A.59) 



6ao J J CCh\K^\ dxdt <fJJ CKx + 54ao j J Qh dxdt 

Qt Qt Qt 

f\Cxx\h\K,\ dxdt <fjj Chi, + 6ao jj CCLh^ dxdt 

Qt Qt Qt 

8ao jj C\Cx\\hx\\Kx\ dxdt < f jj C^^L dxdt + 96ao jj CC^hl dxdt 



Qt 



Qi 



Qt 



Using (A.57)-(A.59) in (A.56) yields 



\H < Y 



/ \ Chi, dxdt + c,jj + C\U + C + ecu {h' + hi) dxdt 

Qt Qt 



Qt 

2\a2\\\w'\ 



(^\Qh dxdt + l\ai\ 1 1 dxdt 



(A.60) 



Qt 



(A.55): 



where ci = max{102ao, 6|ai|} and < e < (v33 — 3)/4. We now consider the term I3 in 
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Integrating by parts, 
/3 + ao Jj^^lx dxdt = f 

Qt 



dxdt + ai ^hl dxdt - ^2 / / {w'C + w"^) h dxdt 



Qt 



Qt 



Qt 



— 2 



J I [l2eC + hi dxdt +Wi\JJ C'hl dxdt 

Qt Qt 

+ 4|a2|||w'||oo JJ ^^^'^'^^^ ^^^^ |a2|||w"||oo JJ d^d^ 

Qt Qt 

Qt 

J J {em+c')hdxdt 



+ 4|a2|(||^i;'||oo + Ik'lloo) 



Qi 



where C2 = max{6ao, |ai|}. 

Using bounds ( |A.60D and ( |A.6lD in ( |A.55D , 



J CGe{h{x,T))dx + ^ jj C^hl^dxdt< j C^Ge{ho,{x)) dx 



+ c-sWeC + C'\U + C + (XL + CWoo + K) dxdt 



Qt 



+ 6|a2|(||^'||oo + |k"||oo)||C'iai+C' 



h dxdt 



Qt 



+ 4|a3|||C'|C:.||oo G,{h)dxdt+\\a^\ // dxdt 



Qt 



Qt 



(A.61) 



(A.62) 



where C3 = max{102ao, 6|ai|} and < £ < (v33 — 3)/4. Taking 5 = in the a priori estimate 
(3.23), using conservation of mass, and assuming < e < minjeo, ("\/33 — 3)/4} where £0 is 
from Lemma 3.1, we deduce from ( |A.62 ) that for all T G [0,Tioc] 



iGe{he{x, T))dx^ / ^Ge{h^e) dx + C 



(A.63) 



n 
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where C > is independent of e > 0. C is determined by Oq, cti, ^2, 03, T^oc, / h^, \Q\ 



w 



w", and on ( and its derivatives. Note that in going from (A. 62) to (A. 63) we dropped the 



ao/2 J J ihl^^^ term because this term is not needed in the rest of the proof. 

We now argue that the e — > limit of the right-hand side of (|A.63[) is finite and bounded 



by allowing us to apply Fatou's lemma to the left-hand side of (A.63), concluding 



^(x) Go{h{x, T)) dx 



h{x,T) 



dx < K < 00 



for every T G [0,T;oc], as desired. (Note that in taking £ ^ we will choose the exact same 
sequence Sk that was used to construct the weak solution h of Theorem [T| Also, in applying 
Fatou's lemma we used the fact that {h = 0} having measure zero in Qti^^ implies {h{-,T)} 
has measure zero in VL.) 

It suffices to show that / ^^^(/loe) J ^Go{ho) < 00 as e — * 0. This uses the Lebesgue 
Dominated Convergence Theorem. First, note that 



1 I _J_ 



22 



Goiz) + 



e 

622 ' 



hence if ho{x) > then 



1 



2{ho(x)+eO) ' 6{ho{x)+eOy 



< 



1 



2hoix) 



+ 



Because has finite entropy (/ Go{hQ) < 00) it is positive almost everywhere in f2. Using 
this and the fact that 6 was chosen so that 6 < 2/5 < 1/2, we have \^{x) Ge{hQ^{x))\ < 
^{x){Go{ho{x)) + c) < G{Go{ho{x)) + c) almost everywhere in x and for all e < Eq. The 
dominating function is in L^, because ho has finite entropy. 

It remains to show pointwise convergence ^(x)G'£(/ioe(x)) C,{x)Go{hQ{x)) almost every- 
where in x: 

\Geihoe{x)) - Go{hoix))\ < \Geihoe{x)) - Goihoe{x))\ 

+ \Go{hoe{x)) - Go{ho{x))\ = ^Tf^ + \Go{hoe{x)) - Go{ho{x))\ 



< 



+ \Goihoeix)) - Goiho{x))\ 



As before, the term ^^/6 goes to zero by the choice of 6*. The term \GQ{hoi;{x)) — Go{hQ{x))\ 
goes to zero for almost every x E fl because Go{z) is continuous everywhere except at z = 0. 

□ 
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B Results used from functional analysis 

Lemma B.l. ([2^) Suppose that X, Y, and Z are Banach spaces, XdF cZ, and X and 
Z are reflexive. Then the embedding {u eLp°{0,T; X) : dtu e L^i(0,T; Z), 1 < pi < oo,i = 
0, 1} (E LPo(0,T; Y) is compact. 

Lemma B.2. ([31^) Suppose that X, Y, and Z are Banach spaces and X (^YcZ. Then the 
embedding {u E L°°(0,T;X) : dtU E Lp{0,T;Z), p > 1} m C{0,T;Y) is compact. 
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